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ABSTRACT 

From the analysis of the near-UV spectrum of the QSO 2206—199, obtained 
with a long series of exposures with STIS on the HST, we deduce a value 
D/H = (1.65 lb 0.35) X 10~^ (Icr error) for the abundance of deuterium in the 
-^abs = 2.0762 damped Lyman a system (DLA) along this sight-line. The velocity 
structure of this absorber is very simple and its neutral hydrogen column density, 
A'^(H I), is accurately known; the error in D/H is mostly due to the limited signal-to- 
noise ratio of the spectrum. Since this is also one of the most metal-poor DLAs, with 
metal abundances ~ 1/200 of solar, the correction due to astration of D is expected to 
be insignificant and the value we deduce should be essentially the primordial abundance 
of deuterium. When all (six) available measurements of D/H in high redshift QSO 
absorbers are considered, we find that the three DLAs — where A^(H I) is measured 
most reliably — give consistently lower values than the three Lyman limit systems. We 
point out that the weighted mean of the DLA measurements, D/H = (2.2 it 0.2) x 10~^, 
yields a baryon density fi^/i^ = 0.025 it 0.001 which is within ~ la of the value deduced 
from the analysis of the CMB angular power spectrum, and is still consistent with the 
present-day D/H and models of Galactic chemical evolution. Future observations of 
D I absorption in other DLAs are needed to establish whether our finding reflects a 
real advantage of DLAs over other classes of QSO absorbers for the measurement of D, 
or is just a statistical fluctuation. 

Subject headings: quasars: absorption lines — quasars individual (Q2206— 199) — 
cosmology: observations — nuclear reactions, nucleosynthesis, abundances 



^Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science 
Institute which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract 
NAS 5-26555. 
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1. Introduction 

Among the light elements created in the Big Bang, deuterium is the one whose abundance 
depends most sensitively on the baryon-to-photon ratio, r], which can in turn be related to the 
cosmological density of baryons Qb, through the known temperature of the cosmic microwave 
background (CMB). It was realized 25 years ago (e.g. Adams 1976) that QSO absorption line 
systems at high rcdshift would provide the most direct avenue to a precise determination of the 
primordial abundance of deuterium, avoiding the need to account for its subsequent destruction 
through successive cycles of star formation in galaxies (astration). However, we have had to wait 
until the advent of 8-10 m telescopes, equipped with efficient high resolution spectrographs, to 
realize those early expectations. 

The reason is simple. The isotope shift between the D and H components of the Lyman series 
amounts to only —82 km s~^; the vast majority of QSO absorption systems span or exceed this 
velocity range so that features due to D (less abundant than H by a factor of a few times 10^) are 

normally blended with other, generally stronger, neutral hydrogen absorption components within 
the same complex. In order to find the rare absorbers with the simplest velocity structure it is 
necessary to reach down the QSO luminosity function to magnitudes which are only accessible with 
the light gathering power of the largest telescopes. Thus, nearly ten years since the commissioning 
of the Keck high resolution echelle spectrograph (HIRES; Vogt 1992) and five years since the first 
detection of D I absorption in a QSO spectrum (Tytler, Fan, & Buries 1996), the available body of 
data on the primordial abundance of D is still frustratingly small, amounting to only five measures 
(Kirkman et al. 2001, also reviewed below). In contrast, determinations of the He/H ratio — far 
more accessible than D/H but much less sensitive to ^Ib — number in the many tens (see Pagel 
2000 for a recent review). 

Until recently, only Lyman limit systems (LLS) — absorbers with neutral hydrogen column 
densities (H I) > 3 x 10^^ cm~^ — had been targeted in searches for interstellar deuterium at high 
redshift. The rarer damped Lyman a systems (DLAs), with A^(H I) > 2 x lO^'^cm"^, potentially 
offer significant advantages over LLS. First and foremost, the neutral hydrogen column density 
can be measured more precisely than is normally the case for LLS. The combination of the width 
of the saturated core and damping wings of the Lyman a absorption line usually constrain A(H I) 
to within 10% or better, even in moderate resolution spectra (e.g. Pettini et al. 1997). Second, 
a knowledge of the degree of metal enrichment of the QSO absorber is required to relate the 
measured D/H ratio to the primordial value. In DLAs the metallicity of the gas can be determined 
straightforwardly and with good precision, without recourse to the uncertain, and often large, 
ionization correction required for LLS (e.g. Viegas 1995). Thirdly, the probability of blending with 
'interloper' H I clouds is much reduced for DLAs, given the steep slope of the A(H I) distribution 
(Peroux et al. 2001), and D absorption should be detectable in many higher order Lyman lines, 
rather than in just the first few transitions in the series. 

In principle, we know of no reason why the velocity structure of DLAs should be systematically 



-3- 



more complex than that of LLS — one could in fact argue that the opposite is more likely to 
be the case, if DLAs are formed in the inner regions of protogalactic disks while LLS arise in 
more extended halo regions (e.g. Steidel 1993). Furthermore, there is mounting evidence, both 
observational (e.g. Pettini et al. 1999 and references therein) and theoretical (e.g. Mo, Mao, &: 
White 1998; Jimenez, Bowen, & Matteucci 1999) that DLAs may trace preferentially galaxies with 
low rates of star formation, where astration of D is also likely to have been low. In any case, we 
believe that DLAs are no less suitable than LLS in the search for high redshift D I absorption; 
they are just rarer and this explains why they have not figured as prominently up to now. 

Among the well studied DLAs, the Zabs = 2.0762 system in the bright (V = 17.3) Zem. = 2.559 
QSO Q2206— 199 is a prime candidate for showing resolved D I absorption. High resolution 
spectroscopy obtained initially with the University College London echelle spectrograph at the 
Anglo-Australian telescope (Pettini &; Hunstead 1990) and later with HIRES on Keck (Prochaska 
Sz Wolfe 1997) showed that in this DLA all the metal lines from neutral gas apparently consist of a 
single, narrow component with a Doppler parameter b = 5.7 km s~^. This is also one of the lowest 
metallicity DLAs known, with different elements being underabundant relative to solar by factors 
between 170 (Si) and 420 (Fe). Given the minimal amount of chemical enrichment experienced by 
the interstellar medium of this absorber, we could reasonably expect its deuterium abundance to 
be very close to the primordial value. 

At the high column densities of DLAs, the D and H components are saturated and blended 
together in the strongest lines of the Lyman series. In order to resolve the two isotopes, it is 
necessary to access high order transitions; at Zabs = 2.0762 all the lines of diagnostic value fall 
below the atmospheric cut-off at 3100 A and are not accessible from the ground. For this reason 
in December 1996 we were awarded 31 orbits of the Hubble Space Telescope to record the near-UV 
spectrum of Q2206— 199 with STIS (the Space Telescope Imaging Spectrograph). For a variety of 
reasons the full set of observations was only completed four years later, in September 2000. In this 
paper we present the resulting co-added spectrum from which we obtain a new measurement of 
the deuterium abundance. In the mean time, D I absorption has been detected in another DLA, 
the Zabs = 3.0249 system in Q0347-393, by D'Odorico, Dessauges-Zavadsky, & Molaro (2001) who 
used the Ultraviolet- Visual Echelle Spectrograph (UVES) on the Very Large Telescope (VLT) of 
the European Southern Observatory. Thus, the total number of D/H measures in high redshift 
QSO absorbers now stands at six; we briefly comment on the dispersion of these values and on 
what it may imply for an improved determination of the primordial abundance of deuterium. 



2. Observations 

The STIS spectra were recorded with the near-UV MAMA detector and G230M grating 
through the 52 x 0.2 arcsec long slit with 31 orbits spread over seven visits between November 13, 
1998 and September 16, 2000. The total integration time was 82 490 s. The grating was set at a 
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central wavelength Ac = 2818 A. With this configuration, the useful wavelength range recorded 
by the near-UV MAMA detector is 2773-2862 A; in the rest frame of the Zabs = 2.0762 DLA this 
corresponds to the interval 901.4-930.4 A which includes all Lyman series lines from Ly7 to the 
Lyman limit. The resolution is 0.19 A (20 km s^^ FWHM) sampled with 2.1 pixels. 

The 31 STIS exposures were processed with the standard CALSTIS pipeline software. We 
extracted one-dimensional spectra from the fluxed and wavelength calibrated frames with the 
IRAF routine apall and then added them together after rebinning to a common wavelength 
scale. We found that the pipeline processing underestimates the scattered light in these low level 
exposures; consequently, we applied an empirical background correction of 5% of the continuum 
level in order to bring the cores of the saturated absorption lines to zero. Figure 1 shows the 
final spectrum, background subtracted, normalized to the QSO continuum, and reduced to the 
rest frame of the Zabs = 2.0762 DLA. At the low light levels of the present observations the 
noise is dominated by the dark count of the NUV MAMA (which turned out to be significantly 
higher than pre-flight expectations and more than one hundred times higher than that of the 
FUV MAMA) and is roughly constant along the spectrum, irrespectively of the QSO signal. Its 
magnitude can be best appreciated from the fluctuations about the zero level in the cores of the 
strong saturated absorption lines near 927.5 and 923 A and below the effective Lyman limit at 
914 A. In the continuum near Ly9 (see below) the signal-to- noise ratio is S/N ~ 10. 



3. Deuterium in the Zabs = 

2.0762 DLA 

As can be appreciated from Figure 1, Q2206— 199 exhibits a crowded spectrum at UV 
wavelengths. The QSO flux is absorbed by a multitude of features due to the Lyman a forest, the 
Zabs = 2.0762 DLA, and a second DLA intercepted in this direction, at Zabs = 1-9205, with much 
stronger and wider interstellar absorption lines (Pettini &: Hunstead 1990; Prochaska & Wolfe 
1997). For the present purposes wc focus on the high order Lyman series of the Zabs = 2.0762 
DLA from Ly7 at 926.2257 A to Lyl3 at 916.429 A; transitions to energy levels with n > 14 are so 
close in wavelength that they become blended with one another eventually producing an effective 
Lyman break near 914.5 A. Fi gure 2 shows on an expanded scale six of the seven Lyman lines 
considered here. Three of them, LylO. Lyll. and Lyl3, are blended with other absorption features, 
as is LyS (shown in Figure 1). The remaining three lines, however, appear to be relatively free of 
contamination and D I absorption is clearly present in Ly7, Ly9, and Lyl2. These are the three 
spectral features on which our measurement of A^(D I) is based. 

All the six Lyman lines in Figure 2 were fitted with Voigt profiles using the VPFIT package.^ 
VPFIT returns the most likely values of redshift z, Doppler width h (km s^^), and column 



^VPFIT is available at http://www.ast.cam.ac.uk/ rfc/vpfit.html 
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density N (cm~^) by minimizing the difference between observed and computed profiles after 
convolution with the appropriate instrumental point spread function (PSF). We used the STIS 
PSF for the G230M grating supplied by the Space Telescope Science Institute; close to the line 
core the function is well approximated by a Gaussian with FWHM = 20 km s^^ but its wings are 
significantly broader. In our fitting procedure we fixed the column density of neutral hydrogen to 
N{}i I) = 2.73 X 10^° cm"^ deduced from the profile of the damped Lyman a line (see §4 below), 
and the number of absorption components to one. All the metal absorption lines from H I gas in 
this DLA are best reproduced by a single absorption component at Zabs = 2.07623 with Doppler 
parameter b = 5.7km s^^ (Prochaska Sz Wolfe 1997; Pettini et al. in preparation). Since the same 
value of b apparently applies to atoms spanning a range of atomic weight from 14 (N) to 56 (Fe), 
this velocity dispersion presumably reflects primarily turbulent, rather than thermal motions. 
Accordingly, we fixed 6turb = 5.0 km and let VPFIT solve for the thermal component, 6th, 
to the H I and D I Lyman lines; the results are collected in Tabic 1. The theoretical absorption 
profiles corresponding to the best solution returned by VPFIT are shown as continuous lines in 
Figure 2. 

We deduce logA'^(D I) = 15.65 it 0.1 (Ifx). The error is dominated by the modest signal-to-noise 
ratio of the STIS spectrum, rather than by systematic uncertainties in the continuum level or the 
details of the fitting procedure. We explored at length the effects of adopting different values 
of 6turb) varying the continuum level and A^(H I) by ±la, assuming a two-component model for 
the absorption lines^, and in all cases VPFIT converged to values of logiV(D I) which were well 
within the ±0.1 dcx uncertainty due to the noise in the spectrum. Figinc 3 illustrates how the 
theoretical profiles obtained by lowering and increasing the value of A^(D I) by 2a compare with 
the observed spectral lines in the two most sensitive cases, Ly9 and Lyl2. The fits shown include 
absorption from adjacent components which are not D I; excluding these components resulted in 
values of A^(D I) which are within the error quoted (although the match to the observed profiles 
was naturally worse). 

We note in passing that b = 14.6 km s~^ (Table 1) corresponds to 6th = 13.7 km for our 
assumed 6turb = 5.0 km s~^ (since 6th 3> 6turb, very similar values of 6th are obtained for all possible 
values of 6turb! from to 5.7km s~^). The implied temperature of the H I gas, T = mb^ /2k where 
m is the atomic mass and k is Boltzmann constant, is T = 11 300 K. This value is essentially the 
same as that deduced in an analogous way by O'Meara et al. (2001) in their determination of 
the abundance of D in the z^bs = 2.536 sub-DLA (7V(H I) = 2.65 x lO^^ cm^^) toward the QSO 
HS 0105-1-1619. Temperatures of ~ 10^ K may appear rather high for H I gas. However, it is 
now well established that high-z DLAs have far higher spin temperatures (deduced by comparing 
21 cm and Lyman a absorption) than H I clouds in the Milky Way. The recent comprehensive 
compilation of values of Tg by Kanekar & Chengalur (2001) includes eight DLAs with lower limits 
ranging from Tg > 800 K to Tg > 4700 K. Both the DLA considered here and the one studied 



^As appropriate to the high ionization gas traced by the Si IV and C IV lines in this DLA 
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by O'Meaxa et al. (2001) have metallicities Z ;$ 1/100^0. In this regime the low cooUng rate 

from metals may well result in interstellar H I temperatures as high as 10^ K. Alternatively, the 
relatively large velocity dispersions deduced by ourselves and O'Meara et al. (2001) may be 
indicative of the presence of additional absorption components, but in this case it would be a 
coincidence that both DLAs yield similar values of T. We stress again, however, that the value 
of 7V(D I) derived above does not depend on the precise details of the velocity structure of the 
absorber because the lines most sensitive to N{Yy I), Ly9 and Lyl2, are both on the linear part of 
the curve of growth. 



4. The Column Density of H I and the Abundance of Deuterium 

Figure 4 shows the normalised spectrum of Q2206— 199 in the region near 3740 A, which 
encompasses the damped Lyman a line at Zabs = 2.07623. This spectrum was obtained by Pettini 
et al. (in preparation) as part of a program to study the abundances of N and O in DLAs. 
The profile of the DLA is best fitted with a column density A/"(H I) = 2.73 x 10^° cm^^ (top 
panel of Figure 4). VPFIT returned a formal error of less than 1%. Even though the strong 
Lyman a line is relatively unblended, we consider this estimate to be over-optimistic. The main 
source of error here is probably the continuum placement within the damped line (from ^ 1200 
to 1235 A — see Figure 4); various trials with different interpolations across this wavelength 
interval showed that ±5 x 10^^ cm~^ is a more realistic assessment of the la error. In any case, 
the error on iV(H I) is clearly much smaller than the error on N(D I) and it is the latter that 
dominates our measure of D/H. For comparison, the lower resolution and S/N spectrum of 
Q2206— 199 obtained by Pettini et al. (1994), with a different detector and telescope, yielded 
A^(H I) = (2.7 ± 0.4) X 10^° cm~^, in excellent agreement with the value deduced here. Evidently, 
the column density of neutral gas in this DLA is known with a satisfactory degree of confidence 
for our present purpose. 

From the above measurements we deduce a value of the abundance of deuterium 
D/H = (1.65 ± 0.35) X 10~^; this is the lowest estimate of this ratio obtained up to now. 



5. Discussion 

At the time when the observations presented here were first proposed, it was still being 
debated whether the primordial abundance of deuterium was a few times 10~^ or one order of 
magnitude higher (Tytler, Fan, & Buries 1996; Songaila et al. 1995). The controversy is now 
largely resolved; with several more low values of D/H reported in the last few years, it seems 
increasingly likely that cases where D/H is apparently greater than 10^'^ are in reality due to 
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contamination by H I at velocities similar to the isotope shift — always a possibility in the Lyman a 
forest (e.g. Kirkman et al. 2001). The finding here of a sixth system with low D/H adds further 
weight to this conclusion. 

In Table 2 we have collected published determinations of D/H in QSO absorption line systems 
at high z; the data are mostly from a similar compilation by O'Mcara et al. (2001 — their Tabic 5) 
augmented by the recent determinations by D'Odorico et al. (2001) and ourselves. The results are 
plotted in Figure 5. Although the set of measurements is still very limited, we comment on two 
points of potential interest. 

All six absorbers are low mctallicity systems, with abundances (as measured by Si) less than 
1/30 of solar ([Si/H] < —1.5). We do not believe that Figure 5o shows a trend of decreasing D/H 
with increasing [Si/H]. The main reason is that in a closed-box model of chemical evolution a 
metallicity Z = —1.5 is reached when, to a rough approximation, only 3% of the gas has been 
processed through stars — and only 3% of the deuterium has been destroyed. Thus an increase 
in metal abundances from —2.5 to —1.5 should result in an imperceptible decrease in the D/H 
ratio through astration. An additional consideration to be aware of is that not all measurements 
of [Si/H] in Figure 5a are equally reliable. While the metallicities of the three DLAs (squares) 
are reasonably secure, the estimates in the three LLS (triangles) rely on much larger ionization 
corrections to account for unobserved ion stages. Thus one should be wary of any apparent 
correlation of D/H with Z when results from LLS and DLAs are mixed together. 

For these reasons it has been assumed until now that in all the cases listed in Table 2 we are 
measuring the primordial abundance of deuterium and that the dispersion between values found in 
different QSO absorbers is just due to observational error. Thus, Tytler et al. (2000) and O'Meara 
et al. (2001) proposed averaging the available determinations to obtain (in the latter, more recent 
work) a most likely value D/H = (3.0 ib 0.4) x 10"^. 

On the other hand, it is interesting that all three DLAs give consistently lower values of 
D/H than the three LLS (Figure 56). In general DLAs allow a more accurate measure of the 
column density of H I than LLS; the shape of the damping wings and the core of the Lyman a 
line normally constrain the allowed values of A^(H I) to a narrow range. This is not always the 
case for LLS where A^(H I) relies more sensitively on the accuracy of the extrapolation of the QSO 
continuum near the Lyman limit and on the details of the velocity structure of the H I gas giving 
rise to saturated Lyman series lines. To illustrate this point, we note that in PKS 1937—1009 the 
initial measure N{li I) = (8.7 ± 1.2) X lO^'^cm-2 by Tytler, Fan, & Buries (1996) was later revised 
by Buries & Tytler (1997) down to N{}1 I) = (7.2 ± 0.3) x lO^^cm'^, while a different analysis 
(Songaila, Wampler, k Cowie 1997) proposed N(R I) < 5 x 10^''' cm'^. In Q1009+299 there is 
the additional problem of partial contamination of the D I absorption by unrelated H I (Buries &; 
Tytler 1998); this is usually less of a concern in the higher column density damped systems. 

We do not intend to be critical of these excellent analyses which have taken great care to 
account for all the relevant factors and arrived at the best estimates allowed by the data. Here 
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we simply speculate, on the basis of the apparent trend in Figure 56, that D/H may have been 
overestimated in the two LLS where deuterium has been detected so far, possibly through an 
underestimate of the column density of neutral hydrogen, and that the most reliable measure of 
the primordial D/H may in fact be provided by the three damped Lyman a systems, where A'^(H I) 
is most secure. Clearly this conjecture will be tested by future observations of additional DLAs; 
nevertheless, it is worthwhile considering briefly its implications, should it turn out to be correct. 

Prom the three DLAs in Table 2 we obtain a weighted mean: 

D/H = (2.2 ±0.2) X 10"^ (1) 

where the weights used are inversely proportional to the square of the error on each measurement 
and 0.2 x 10~^ is the la error on the weighted mean. Using the analytic expression by 
Buries, Nollett, & Turner (2001), the primordial abundance of deuterium in eq. (1) implies a 
baryon-to-photon ratio: 

T] = (6.8 ± 0.4) X 10"^° (2) 

which in turn yields a baryon density 

Ob/i^ = 0.025 ±0.001 (3) 

using the scaling of fis/i^ with rj given by Buries et al. (2001). 

We note that a primordial abundance of deuterium as low as D/H = (2.2 ± 0.2) x 10~^ 
can still be reconciled with the present-day value (D/H)jgj^ = (1.6 ± 0.1) x 10~^ in the Galactic 
interstellar medium (Tytler et al. 2000 and references therein), particularly when the unrealistic 
(for the Milky Way at least) assumption of closed-box chemical evolution is relaxed to include 
outflows and infall of unprocessed gas (e.g. Edmunds 1994; Tosi et al. 1998). It is also possible 
that (D/H)jgj^,j is less than 1 x 10~^, given two recent lower determinations (Vidal-Madjar et al. 
1998; Jenkins et al. 1999). Turning to other light elements produced in the Big Bang, the value 
of rj in eq. (2) implies a primordial ^He abundance Yp = 0.249 which is still within the range 
of current estimates (Pagel 2000), but a high value of the primordial Li/H= 6.0 x 10~^° which 
exacerbates the well-known lithium problem (Ryan et al. 2000; Buries et al. 2001). 

In the last few months there has been considerable discussion of a possible conflict between 
the values of JIb deduced from the abundances of the light elements interpreted within the 
framework of Big-Bang nucleosynthesis on the one hand, and from the analysis of the angular 
power spectrum of the CMB on the other. Recent measurements of the latter have been analyzed 
to give a best fitting 

neh^ ^ 0.0S2t^^Zl (4) 

where the error is again la (Jaffe et al. 2001). The 2a difference from the baryon density implied 
by existing measurements of D/H as summarized by Tytler et al. (2000) and O'Meara et al. 
(2001) has prompted (not for the first time) speculations that 'new physics' — or at least previously 
unaccounted for processes (e.g. Naselsky et al. 2001) — may be required to reconcile these two 
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determinations. On the other hand, if the low values of D/H found in damped Lyman a systems 
are indeed representative of the true primordial abundance of D, as speculated here, the apparent 
conflict is reduced considerably. The two estimates in eqs. (3) and (4) are nearly within la of each 
other and the history of astronomical measurements teaches us that la errors almost invariably 
turn out to be over-optimistic in the light of subsequent improvements in the data. 

The key question is whether the systematically low values of D /H in the three DL As in Figure 
5 are a significant result or just a statistical accident. This question can only be answered by more 
observations of DLAs with the characteristics required to detect and measure D I absorption. 
Fortunately the prospects for finding such valuable test cases are excellent, given the spectacular 
increase in the numbers of known QSOs brought about by large-scale searches such as the Sloan 
digital sky survey and 2dF QSO survey. A significant investment in observing time will be required 
to follow up these QSOs with higher resolution spectroscopy. However, it is now clear that only 
by assembling a moderately large sample of D /H measurements at high redshift will it be possible 
to differentiate between systematic and random errors and eventually determine the primordial 
abundance of deuterium. 



6. Epilogue 

Since this paper was submitted for publication there have some important developments 
of relevance to the discussion above. Three new estimates oi have been announced, all 
based on observations of the CMB power spectrum. From the analysis of the full dataset of the 
BOOMERANG experiment, Netterfield et al. (2001) and de Bernardis et al. (2001) deduced 

= 0.0211q;qq| (la errors). The first results from the Degree Angular Scale Interferometer 
(DASI) also give 9.Bh? = 0.022to;oo3 (Pryke et al. 2001). A new treatment of the MAXIMA-I 
measurements favours a value which is higher, fi^/i^ = 0.0325 it 0.007, but still consistent with 
the other two within the errors (Stompor et al. 2001). Thus one could reasonably argue that the 
'problem' has been solved and that fi^ is now known to within ^ 10%, since the weighted average 
of all five detections in Table 2, D/H = (2.6 ± 0.2) x 10"^, implies J^b/i^ = 0.022 ± 0.001, in 
excellent agreement with the latest CMB mesurements. It just remains to be seen how this newly 
found concordance stands the test of future measurements. 

A second recent development is the re-analysis by Levshakov et al. (2001) of the absorption 
lines in the Zabs = 3.0249 DLA in Q0347-383 previously studied by D'Odorico et al. (2001). 
By including in the complex velocity structure of this DLA a narrow (b = 3km s~^) component 
whose presence is suggested by S II and H2 absorption lines, Levshakov et al. propose that the 
abundance of D should be revised upwards from (D/H) = (2.25 it 0.65) xlO~^ (the value derived 
by D'Odorico et al. and quoted in Table 2) to (D/H) = (3.2 ± 0.4) xlO"^. If this is indeed the 
case, then the evidence for a systematic difference in the values of (D/H) between Lyman limit 
and damped Lyman a systems is made significantly weaker. As emphasized above, only further 
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searches for D in absorption systems with a simple velocity structure, such as that presented here, 
will clarify this point. 

In any case, as far as primordial nucleosynthesis is concerned, one of the most urgent goals 
now is identifying the origin of the dispersion in the values of D/H in Table 2 and Figure 5. The 
simplest explanation is that the errors of the individual measurements have been underestimated. 
However, it is worrying that a similarly wide dispersion in D/H can be found in the local ISM, 
where the factor of ~ 2 variation indicated by earlier work (Sonneborn et al. 2000 and references 
therein) is apparently confirmed by the latest, high precision, observations with the Far Ultraviolet 
Spectroscopic Explorer (Moos et al. 2001). At both high and low redshift, we do not have a 
satisfactory explanation for differing D/H ratios in gas of similar chemical composition. Resolving 
this puzzle is important for two reasons. First, it dents our confidence in the general framework of 
Big-Bang nucleosynthesis. Second, it thwarts future attempts to use the abundance of D, together 
with those of elements which are manufactured — rather than destroyed — by stars, to track the 
chemical evolution of galaxies over the age of the universe. Such attempts seem somewhat futile 
at present given that the total amount of astration of D over the Hubble time is comparable to 
the scatter in the determination of D/H at any redshift. 

We arc most grateful to Bob Carswell, who kindly performed for us the line profile fits 
presented in this paper, and to Bernard Pagel for helpful comments. MP would like to acknowledge 
J. Bergeron, P. Petitjean, S. Ellison, and C. Ledoux, his collaborators in the VLT program to 
study DLAs; the data shown in Figure 4 were obtained as part of that study. 
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Table 1. BEST FITTING PARAMETERS FOR H I and D I ABSORPTION AND la ERRORS 



Ion 


N (cm"2) 


b (km s ^) 


z 


H I 


(2.73 ± 0.05) X lO^O" 


14.6 ± 0.2^ 


2.076234 ± 0.000002 


D I 


(4.5 ± 1) X 10^5 


10.6^ 


2.076234'^ 



^Determined from the profile of the damped Lyman a line 

'^The Dopplcr parameter b has both thermal and turbulent components (6^ = b"^^ + ?'turb)' ^turb was 
fixed at 5 km (see text) 

'^Fixed to be 6(D I)^ = b^^^^ + 1/2 x 6th (H I)^ 
'^Fixed to be the same as ^;(H I) 



Table 2. SUMMARY OF D/H MEASUREMENTS AT HIGH REDSHIFT 



QSO 




iV(H I) (cm-2) 


[Si/H]" 


(D/H)±lcT (10^5) 


Ref. 


Q0130-403 


2.799 


4.6 X 10^^ 


-2.6 


< 6.8 


O'Meara et al. (2001) 


Q1009+299 


2.504 


2.5 X 10^^ 


-2.53 


4.0 ±0.65 


O'Meara et al. (2001) 


PKS 1937-1009 


3.572 


7.2 X 10^^ 


-2.26 


3.25 ±0.3 


O'Meara et al. (2001) 


HS 0105+1619 


2.536 


2.6 X 10^9 


-2.0^ 


2.5 ± 0.25 


O'Meara et al. (2001) 


Q2206-199 


2.0762 


2.73 X 10^0 


-2.23^ 


1.65 ± 0.35 


This paper 


Q0347-383 


3.025 


4.3 X 10^° 


-1.53'^ 


2.25 ± 0.65 


D'Odorico et al.(2001) 



'^In the usual notation, [Si/H] = log(Si/H) — log(Si/H)Q. We chose Si as an indicator of the overall 
mctallicity Z because this is the element for which most measurements are available for the absorption 
systems considered here 

'^This is actually [0/H]. The upper limit [Si/H] < -1.85 determined by O'Meara et al. (2001) is 
consistent with the oxygen abundance 

'^Prochaska & Wolfe (1997) 
•^Prochaska k Wolfe (1999) 
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Fig. 1.— Near-UV (2805.5-2857.8 A) STIS spectrum of Q2206-199 normalised to the QSO 
continuum and reduced to the rest frame of the z^bs = 2.0762 DLA. The spectral resolution is 
20 km s"^ FWHM and S/N ~ 10; the total exposure time was 82 490 s. The locations of high 
order Lyman lines of H I and D I, from Ly7 to Lyl3, are indicated by vertical tick marks above 
the spectrum. 
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Fig. 2. — Residual intensity vs. rest wavelength (A) for high order Lyman lines in the 2;abs — 

2.0762 

DLA. The histograms are the data, while the thin continuous lines show the theoretical absorption 
profiles for H I and D I corresponding to the best fit solution returned by VPFIT (see Table 1). 
The important lines here are Ly7, Ly9, and Lyl2 which are sufficiently free of contamination to 
allow a determination of the column density of D I. The fits shown for these three lines include 
adjacent, unrelated, absorption features. For LylO, Lyll, and Lyl3 we show only the contributions 
of the H I and D I components to the complex blends present near their wavelengths. Although 
there is absorption at the expected position of D I in each of these lines, the blends in these three 
cases are dominated by contaminating features, as evidenced by the fact that the optical depths do 
not match the relative /-values of LylO, Lyll, and Lyl3. 
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Fig. 3. — Comparison between observed (histogram) and computed (continuous line) profiles for 
logA'^(D I) = 15.65 lb 0.2. This range corresponds to twice the la error on A^(D I) returned by 
VPFIT. 
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Fig. 4.— Damped Lyman a line in the Zabs = 2.0762 DLA in Q2206-199 obtained with UVES on 
the VLT (reproduced from Pettini et al. in preparation). The resolution is 7km s~^ FWHM. The 
fit shown in the top panel is for A^(H I) = 2.73 x lO^'^ cm~^. The lower panel shows the damped line 
on an expanded scale together with three fits, corresponding to A^(H I) = (2.73 it 0.5) x 10^'^ cm~-^. 
The range shown corresponds to 10 times the estimated error on A^(H I). 
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Fig. 5. — Available measurements of the abundance of deuterium in QSO absorption systems at 
high redshift from Table 2 — triangles: Lyman limit systems; circles: damped Lyman a systems. 
The top panel shows D/H as a function of metallicity, as measured by the [Si/H] ratio. In the lower 
panel the deuterium abundance is plotted against the column density of neutral gas. Damped 
Lyman a systems seem to yield systematically lower values of D/H than Lyman limit systems. 



